library(readxl)
library(tidyverse)
sep07 <- read.csv("../data/clean_data/npoc_7sep21.csv")%>%
mutate(Day = as.Date("2021-09-07"))
sep08 <- read.csv("../data/clean_data/npoc_8sep21.csv")%>%
mutate(Day = as.Date("2021-09-08"))
sep09 <- read.csv("../data/clean_data/npoc_9sep21.csv")%>%
mutate(Day = as.Date("2021-09-09"))
sep14 <- read.csv("../data/clean_data/npoc_14sep21.csv")%>%
mutate(Day = as.Date("2021-09-14"))
sep15 <- read.csv("../data/clean_data/npoc_15sep21.csv")%>%
mutate(Day = as.Date("2021-09-15"))
sep16 <- read.csv("../data/clean_data/npoc_16sep21.csv")%>%
mutate(Day = as.Date("2021-09-16"))
sep17 <- read.csv("../data/clean_data/npoc_17sep21.csv")%>%
mutate(Day = as.Date("2021-09-17"))
sep20 <- read.csv("../data/clean_data/npoc_20sep21.csv")%>%
mutate(Day = as.Date("2021-09-20"))
sep21 <- read.csv("../data/clean_data/npoc_21sep21.csv")%>%
mutate(Day = as.Date("2021-09-21"))
full = sep07%>%
full_join(sep08)%>%
full_join(sep09)%>%
full_join(sep14)%>%
full_join(sep15)%>%
full_join(sep16)%>%
full_join(sep17)%>%
full_join(sep20)%>%
full_join(sep21)%>%
filter(!grepl("spk|dup",Sample))
Joining, by = c("Sample", "NPOC", "TNb", "Day")
Joining, by = c("Sample", "NPOC", "TNb", "Day")
Joining, by = c("Sample", "NPOC", "TNb", "Day")
Joining, by = c("Sample", "NPOC", "TNb", "Day")
Joining, by = c("Sample", "NPOC", "TNb", "Day")
Joining, by = c("Sample", "NPOC", "TNb", "Day")
Joining, by = c("Sample", "NPOC", "TNb", "Day")
Joining, by = c("Sample", "NPOC", "TNb", "Day")
full$Sample_day <- as.numeric(sub("[a-z]+[0-9]+d","",sub("([a-z](_.+)*$)","",full$Sample)))
full$Rep <- as.numeric(substr(full$Sample,4,4))
full$filt <- sub("[a-z]+[0-9]+[a-z]+[0-9]+","",sub("_.+","",full$Sample))
full$Treatment <- substr(full$Sample,1,2)
full_conc = full
full_conc$Treatment[full_conc$Treatment=="aa"]<-"Anoxic"
full_conc$Treatment[full_conc$Treatment=="oo"]<-"Oxic"
full_conc$Treatment[full_conc$Treatment=="oa"]<-"Oxic to anoxic"
full_conc$Treatment[full_conc$Treatment=="ao"]<-"Anoxic to oxic"
fer = full_conc
fer$filt[fer$filt=="s"]<-"Filtered"
fer$filt[fer$filt=="t"]<-"Unfiltered"
fer$Status <- "Oxic"
fer$Status[fer$Treatment %in%c("Anoxic to oxic","Anoxic")]<- "Anoxic"
write.csv(fer,"../data/processed_data/TOC_formatted.csv", row.names = F)
full_conc$Treatment[full_conc$Treatment=="Anoxic to oxic"&full_conc$Sample_day<=13]<-"Anoxic"
full_conc$Treatment[full_conc$Treatment=="Oxic to anoxic"&full_conc$Sample_day<=13]<-"Oxic"
f = fer%>%
filter(!is.na(Treatment))%>%
ggplot(aes(x = Sample_day, y = NPOC, color = Treatment))+
geom_point()+
ylab("NPOC (mg/L)")+
geom_vline(xintercept = 13.5)+
scale_color_manual(values = c("#175676","#D58936","#A44200","#4BA3C3")) +
facet_grid(Status~filt)+
theme_bw()
fer%>%
filter(!is.na(Treatment))%>%
ggplot(aes(x = Sample_day, y = TNb, color = Treatment))+
geom_point()+
ylab("Total bound nitrogen (mg/L)")+
xlab("Day of experiment")+
geom_vline(xintercept = 13.5)+
scale_color_manual(values = c("#175676","#D58936","#A44200","#4BA3C3")) +
facet_grid(Status~filt)+
theme_bw()

library(plotly)
ggplotly(f)
fer%>%
filter(Sample_day == 9,
filt == "Filtered",
Treatment == "Oxic")
fer
f

LS0tCnRpdGxlOiAiUiBOb3RlYm9vayIKb3V0cHV0OiBodG1sX25vdGVib29rCi0tLQoKYGBge3J9CmxpYnJhcnkocmVhZHhsKQpsaWJyYXJ5KHRpZHl2ZXJzZSkKCnNlcDA3IDwtIHJlYWQuY3N2KCIuLi9kYXRhL2NsZWFuX2RhdGEvbnBvY183c2VwMjEuY3N2IiklPiUKICBtdXRhdGUoRGF5ID0gYXMuRGF0ZSgiMjAyMS0wOS0wNyIpKQpzZXAwOCA8LSByZWFkLmNzdigiLi4vZGF0YS9jbGVhbl9kYXRhL25wb2NfOHNlcDIxLmNzdiIpJT4lCiAgbXV0YXRlKERheSA9IGFzLkRhdGUoIjIwMjEtMDktMDgiKSkKc2VwMDkgPC0gcmVhZC5jc3YoIi4uL2RhdGEvY2xlYW5fZGF0YS9ucG9jXzlzZXAyMS5jc3YiKSU+JQogIG11dGF0ZShEYXkgPSBhcy5EYXRlKCIyMDIxLTA5LTA5IikpCnNlcDE0IDwtIHJlYWQuY3N2KCIuLi9kYXRhL2NsZWFuX2RhdGEvbnBvY18xNHNlcDIxLmNzdiIpJT4lCiAgbXV0YXRlKERheSA9IGFzLkRhdGUoIjIwMjEtMDktMTQiKSkKc2VwMTUgPC0gcmVhZC5jc3YoIi4uL2RhdGEvY2xlYW5fZGF0YS9ucG9jXzE1c2VwMjEuY3N2IiklPiUKICBtdXRhdGUoRGF5ID0gYXMuRGF0ZSgiMjAyMS0wOS0xNSIpKQpzZXAxNiA8LSByZWFkLmNzdigiLi4vZGF0YS9jbGVhbl9kYXRhL25wb2NfMTZzZXAyMS5jc3YiKSU+JQogIG11dGF0ZShEYXkgPSBhcy5EYXRlKCIyMDIxLTA5LTE2IikpCnNlcDE3IDwtIHJlYWQuY3N2KCIuLi9kYXRhL2NsZWFuX2RhdGEvbnBvY18xN3NlcDIxLmNzdiIpJT4lCiAgbXV0YXRlKERheSA9IGFzLkRhdGUoIjIwMjEtMDktMTciKSkKc2VwMjAgPC0gcmVhZC5jc3YoIi4uL2RhdGEvY2xlYW5fZGF0YS9ucG9jXzIwc2VwMjEuY3N2IiklPiUKICBtdXRhdGUoRGF5ID0gYXMuRGF0ZSgiMjAyMS0wOS0yMCIpKQpzZXAyMSA8LSByZWFkLmNzdigiLi4vZGF0YS9jbGVhbl9kYXRhL25wb2NfMjFzZXAyMS5jc3YiKSU+JQogIG11dGF0ZShEYXkgPSBhcy5EYXRlKCIyMDIxLTA5LTIxIikpCgpmdWxsID0gc2VwMDclPiUKICBmdWxsX2pvaW4oc2VwMDgpJT4lCiAgZnVsbF9qb2luKHNlcDA5KSU+JQogIGZ1bGxfam9pbihzZXAxNCklPiUKICBmdWxsX2pvaW4oc2VwMTUpJT4lCiAgZnVsbF9qb2luKHNlcDE2KSU+JQogIGZ1bGxfam9pbihzZXAxNyklPiUKICBmdWxsX2pvaW4oc2VwMjApJT4lCiAgZnVsbF9qb2luKHNlcDIxKSU+JQogIGZpbHRlcighZ3JlcGwoInNwa3xkdXAiLFNhbXBsZSkpCgpmdWxsJFNhbXBsZV9kYXkgPC0gYXMubnVtZXJpYyhzdWIoIlthLXpdK1swLTldK2QiLCIiLHN1YigiKFthLXpdKF8uKykqJCkiLCIiLGZ1bGwkU2FtcGxlKSkpCmZ1bGwkUmVwIDwtIGFzLm51bWVyaWMoc3Vic3RyKGZ1bGwkU2FtcGxlLDQsNCkpCmZ1bGwkZmlsdCA8LSBzdWIoIlthLXpdK1swLTldK1thLXpdK1swLTldKyIsIiIsc3ViKCJfLisiLCIiLGZ1bGwkU2FtcGxlKSkKZnVsbCRUcmVhdG1lbnQgPC0gc3Vic3RyKGZ1bGwkU2FtcGxlLDEsMikKZnVsbF9jb25jID0gZnVsbAoKZnVsbF9jb25jJFRyZWF0bWVudFtmdWxsX2NvbmMkVHJlYXRtZW50PT0iYWEiXTwtIkFub3hpYyIKZnVsbF9jb25jJFRyZWF0bWVudFtmdWxsX2NvbmMkVHJlYXRtZW50PT0ib28iXTwtIk94aWMiCmZ1bGxfY29uYyRUcmVhdG1lbnRbZnVsbF9jb25jJFRyZWF0bWVudD09Im9hIl08LSJPeGljIHRvIGFub3hpYyIKZnVsbF9jb25jJFRyZWF0bWVudFtmdWxsX2NvbmMkVHJlYXRtZW50PT0iYW8iXTwtIkFub3hpYyB0byBveGljIgoKZmVyID0gZnVsbF9jb25jCmZlciRmaWx0W2ZlciRmaWx0PT0icyJdPC0iRmlsdGVyZWQiCmZlciRmaWx0W2ZlciRmaWx0PT0idCJdPC0iVW5maWx0ZXJlZCIKCmZlciRTdGF0dXMgPC0gIk94aWMiCmZlciRTdGF0dXNbZmVyJFRyZWF0bWVudCAlaW4lYygiQW5veGljIHRvIG94aWMiLCJBbm94aWMiKV08LSAiQW5veGljIgoKd3JpdGUuY3N2KGZlciwiLi4vZGF0YS9wcm9jZXNzZWRfZGF0YS9UT0NfZm9ybWF0dGVkLmNzdiIsIHJvdy5uYW1lcyA9IEYpCgpmdWxsX2NvbmMkVHJlYXRtZW50W2Z1bGxfY29uYyRUcmVhdG1lbnQ9PSJBbm94aWMgdG8gb3hpYyImZnVsbF9jb25jJFNhbXBsZV9kYXk8PTEzXTwtIkFub3hpYyIKZnVsbF9jb25jJFRyZWF0bWVudFtmdWxsX2NvbmMkVHJlYXRtZW50PT0iT3hpYyB0byBhbm94aWMiJmZ1bGxfY29uYyRTYW1wbGVfZGF5PD0xM108LSJPeGljIgoKZiA9IGZlciU+JQogIGZpbHRlcighaXMubmEoVHJlYXRtZW50KSklPiUKICBnZ3Bsb3QoYWVzKHggPSBTYW1wbGVfZGF5LCB5ID0gTlBPQywgY29sb3IgPSBUcmVhdG1lbnQpKSsKICBnZW9tX3BvaW50KCkrCiAgeWxhYigiTlBPQyAobWcvTCkiKSsKICBnZW9tX3ZsaW5lKHhpbnRlcmNlcHQgPSAxMy41KSsKICBzY2FsZV9jb2xvcl9tYW51YWwodmFsdWVzID0gYygiIzE3NTY3NiIsIiNENTg5MzYiLCIjQTQ0MjAwIiwiIzRCQTNDMyIpKSArCiAgZmFjZXRfZ3JpZChTdGF0dXN+ZmlsdCkrCiAgdGhlbWVfYncoKQoKZmVyJT4lCiAgZmlsdGVyKCFpcy5uYShUcmVhdG1lbnQpKSU+JQogIGdncGxvdChhZXMoeCA9IFNhbXBsZV9kYXksIHkgPSBUTmIsIGNvbG9yID0gVHJlYXRtZW50KSkrCiAgZ2VvbV9wb2ludCgpKwogIHlsYWIoIlRvdGFsIGJvdW5kIG5pdHJvZ2VuIChtZy9MKSIpKwogIHhsYWIoIkRheSBvZiBleHBlcmltZW50IikrCiAgZ2VvbV92bGluZSh4aW50ZXJjZXB0ID0gMTMuNSkrCiAgc2NhbGVfY29sb3JfbWFudWFsKHZhbHVlcyA9IGMoIiMxNzU2NzYiLCIjRDU4OTM2IiwiI0E0NDIwMCIsIiM0QkEzQzMiKSkgKwogIGZhY2V0X2dyaWQoU3RhdHVzfmZpbHQpKwogIHRoZW1lX2J3KCkKCmxpYnJhcnkocGxvdGx5KQpnZ3Bsb3RseShmKQpmZXIlPiUKICBmaWx0ZXIoU2FtcGxlX2RheSA9PSA5LAogICAgICAgICBmaWx0ID09ICJGaWx0ZXJlZCIsCiAgICAgICAgIFRyZWF0bWVudCA9PSAiT3hpYyIpCgpmZXIKZgpgYGAKCg==